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Complementary nature of coercivity enhancement and exchange bias is generalized from 
the layered systems to ferro-antiferromagnet (F-AF) exchange coupled systems with 
arbitrary configurations and is proved based on the coherent rotation of the F 
magnetization. In the proof, the effect of F-AF coupling is absorbed into the anisotropy 
of the F part, resulting in an arbitrary anisotropy for the F part. The proof starts from a 
general discussion on the initial susceptibility of a single domain particle. Then the 
fundamental correlation between the maximal initial susceptibility and the critical field 
along an arbitrary easy direction of a single domain particle, at which the magnetization 
becomes unstable, is discussed. Finally, the difference of the initial switching field and 
the actual switching field along the arbitrarily chosen easy direction is discussed. 
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I. Introduction 



Exchange anisotropy was discovered by Meiklejohn and Bean in 1956. It originates 
from the exchange interactions between the spins of the ferromagnet (F) and the spins of 
the antiferromagnet (AF) at their interface. It is a complicated phenomenon that depends 
on a variety of factors including the ferro-antiferromagnet (F-AF) interface nature and the 
properties of the antiferromagnet. The difficulty of monitoring the F-AF interface and 
the lack of knowledge about the properties of antiferromagnets make the study of this 
anisotropy rather difficult. Furthermore, because of the F-AF coupling, the AF spins are 
not rigid during F magnetization processes. As a result, the effective anisotropy of the F 
part due to the F-AF coupling is not easy to predict. In fact, study of the NigoFe2o/FeMn 
F-AF bilayers showed that the angular dependences of the exchange bias field and the 
coercivity of the Ni 80 Fe2o were sensitive to the film growth direction 2 and were different 
for epitaxial films and textured films . 

Despite these difficulties, there are some general features for the F-AF systems that are 
independent of the details of the systems. For example, it was found that when a 
ferromagnet is coupled to an antiferromagnet through their interface, its magnetization 
behavior is greatly altered; namely, the magnetization tends to be more uniform upon 
reversal compared with the uncoupled single ferromagnet, a law first found by Fujiwara 
et al 4 . This fact is evidenced by the experimental results of the F-AF exchange bilayers 
from the AMR measurements 4 , MFM observations 5 , torque curves 6 , angular dependences 



of the exchange bias field and the coercivity of the pinned layer , etc. But the most 
important experimental result showing the coherent rotation of the F magnetization is the 
so-called complementary nature of coercivity enhancement and exchange bias 4 8 9 10 , 
which was proposed by Fujiwara 4 . For an F-AF exchange bilayer, one can define the 
intrinsic pinning field H p0 by M/xt„±, where M s is the saturation magnetization of the F 
layer and XirA- is the initial susceptibility transverse to the stable direction of the F 
magnetization in the absence of external field. This definition is a natural generalization 
of the anisotropy field for a single F layer. It was found that H p q is approximately equal 
to the switching field of the F magnetization along the stable direction of the F 
magnetization in the absence of external field, denoted by H sw ~, which is the sum of the 
total coercivity H c of the F layer and the exchange bias field H eh . This reveals that in the 
F-AF bilayers both the coercivity enhancement of the F layer and the exchange bias 
originate from the F-AF coupling, and that the magnetization reversal process of the F 
layer is dominated by coherent rotation. Theoretically it was shown that , in magnitude, 
H p o is precisely equal to the F layer's initial switching field along the referred stable 
direction H sw q, defined as the reverse field at which the F magnetization starts to be 
unstable. The difference between the actual switching field and the initial switching field 
will be discussed in more detail in section IV. To avoid confusion, from now on, we use 
the term of the complementary nature of coercivity enhancement and exchange bias in 
the F-AF bilayers to refer to the exact equivalence relationship between H p o and H sw q. 
The original mathematical proof assumed the magnetization to rotate in plane 8 . 
However, the key to our proof was only the coherence of the F magnetization. 



Exchange anisotropy was first found in a compact of fine particles of cobalt with a 
cobaltous oxide shell 1 , for which the F-AF coupling occurs at the Co-CoO interface. By 
nature it was not a layered system. Although the complementary nature was found in F- 
AF bilayers, one would expect a similar relationship in a general F-AF system with 
arbitrary configuration, such as the Co-CoO system and a trilayer AF-F-AF system. 
However, the generality of spin configuration makes it necessary to redefine H p0 . For a 
general F-AF system, we may define H p o as M/xin.max, where Xin,max is the maximum initial 
susceptibility (over all the directions in the 3-d space) attained in some direction 
perpendicular to the stable direction of the magnetization in the absence of external field. 
It is the purpose of this paper to prove, based on the coherence of the F magnetization, 
H p0 = H swQ _ for a general F-AF system, with the definition of H sw0 ' being the same as for 

the layered systems. The strategy for the proof of this particular result is that the effect of 
the F-AF coupling can be absorbed into the anisotropy energy of the F part. That is, the 
overall anisotropy of the F part includes its intrinsic anisotropy as well as the anisotropy 
induced by the F-AF coupling. By doing that the complication of microscopic F-AF 
interactions at the interface can be technically avoided and the generality of the 
discussion is achieved. In fact the same approach allows us to discuss any type of F-AF 
system. Note that the F-AF coupling induced anisotropy is not limited to a unidirectional 
anisotropy plus a uniaxial anisotropy 2 . In order for our discussion to be valid for any 
ferromagnetic materials for the F part of system, we do not assume any specific form for 
the intrinsic anisotropy of the F part, either. With these factors being taken into account, 
the overall anisotropy of the F part in our proof will be of arbitrary form. The proof for 



the complementary nature of the coercivity enhancement and exchange bias in F-AF 
coupled systems thus becomes a problem of an uncoupled single domain particle with 
arbitrary anisotropy. 

For a general 3-d F-AF system for which the initial state is given, the magnetic state 
includes the AF domain structure and is path dependent due to the switching of the AF 
grains during the F magnetization processes. The overall anisotropy of the F part does 
not exist for the whole state space. However, in the initial susceptibility measurement the 
F magnetization only undergoes small deviation from the initial stable direction. Thus if 
the physical processes for perturbations of the F magnetization around the initial stable 
direction are reversible, the overall anisotropy is well defined in the neighborhood of the 
initial stable direction. For a 2-d F-AF system, that constraint can be relaxed because the 
F magnetization rotates in plane. During its switching process it has only one path to 
undergo and always rotates in one sense. Therefore for a 2-d F-AF system, although 
hysteresis occurs, the overall anisotropy is well defined except at switching of the F 
magnetization, and it is thought to be continuous during the continuous rotation of the F 
magnetization. Thus the overall anisotropy is well defined in the initial susceptibility 
measurement for a 2-d system. For the following discussions to be valid, we require that 
the second derivatives of the overall anisotropy exist in a neighborhood of the initial F 
magnetization direction. 



In the following we start by a general discussion of the initial susceptibility for a single 
domain particle in section II. In section III we show H p0 =H . In section IV we 

discuss the difference between the initial switching field and the actual switching field. 
We conclude in section V. 



II. Initial susceptibility 



For any single domain particle, the direction of the magnetization M determines the 
magnetic state of the particle. Without loss of generality, we assume M initially points to 
the positive x direction, one of the easy directions of the particle. In an F-AF system, that 
corresponds to one of the stable directions in the absence of external field, which could 
be different from the pinned direction (see section IV). At this stage, the other two 
coordinate axes can be specified arbitrarily to form a standard 3-d rectangular coordinate 
system. Suppose that an external field H is then applied and the magnetization rotates to 
a new direction. Spherical coordinates (6, <p) and (a, /?) are used to specify the directions 
of H and M, respectively. Fig. 1 sketches the directions of H and M. 
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Let E(a,P) be the anisotropy energy density of the particle. Denote dE/da and d E/da , 
the first two partial derivatives of E(afi) with respect to a, by E a and E aa , etc. The 
assumption for an easy direction at the positive x direction (i.e., a = K I '2, f5 = 0) ensures 
the existence of the partial derivatives of E(a,P) at a = 7t I '2, f3 = . Furthermore, we 



have E a (;z72,0) = 0, E p {nl2,Q) = 0, E aa (x/2,0)>0, E^/2,0) > and 
£ m (7ll2$)E pp (tt/2,0) > [E a/3 (n 7 2,0)] 2 . For simplicity, we can properly select y-axis 
and z-axis so that E a/3 (xl2,0) = (the proof for that is in appendix A). 

The free energy density F of the particle includes the anisotropy energy and Zeeman 
energy 

F(a,(3) = E(a,/3)-M s H[sm0smacos(/3-<p) + cos0cosa], (1) 
where M s is the saturation magnetization. 

At any stage of initial susceptibility measurement, the first derivatives of the free energy 
with respect to a and (3 vanish, yielding 

E a -M s H[sm0cosacos(/3-(p)-cos0sma] = O, (2) 

E /j +M s Hsm0smasm(j3-p) = O. (3) 

To calculate susceptibility, we take a variation $H of the external field along its original 
direction. This results in a change of the magnetization direction, denoted by da and <S/?. 
Thus the new direction of the magnetization is a = da + 71 1 2, /? = 8f3 . From the 
variations of (2) and (3), setting a = 7tl2, f3 = and H - for initial susceptibility, we 
can solve the resulting variational equations for Sa/SH and Sp/SH to obtain 



8al8H = -M s cosOIE aa {nl2,0), 
8/3 1 8H = M s sin 6 sin (p I E pp (n 1 2,0) . 



From the definition of the initial susceptibility Xin we obtain 

%in =8M n l8H 

= M s [- cos 08a I SH + sin sin <pSj3 1 SH ] 

_ M s 2 [E aa {n 1 2,0) sin 2 9 sin 2 p + £^ <>/ 2,0) cos 2 6>] 

M 2 [P (^/ 2,0)e 2 „, + 2,0)e 2 „ z ] 



(4) 



^(^r/2,0)^(^/2,0) 
where e Hy - sin # sin q> and e Hz = cos# are the direction cosines of the external field H, 
and 

M h = M s [sin 6>sin acos(/3 -<p) + cos 6 cos a] 
is the parallel component of the magnetization along H. 



From (4), we see that/ in has a great deal of symmetries. Namely, it is symmetric with 
respect to the x-y, y-z, and x-z planes. This symmetry can be easily verified from the 
definitions of eHy and em- Furthermore, Xin depends on E(a,[}) only through E aa (nl2,Q>) 
and E m (x/2,0). 



The above properties of are reasonable since Xin involves the response of the 
magnetization to the perturbation of an infinitesimal field. Thus it depends only on the 
local anisotropy as the magnetic state is being perturbed. The symmetry of the local 
anisotropy determines the symmetry of Xin- hi our case, the perturbed state is around 
a - 71 12, /? = . Near this direction, the anisotropy energy density £ at a = 7tl2 + da , 
f3 = 8f5 is (to the second order) 

E(a,jB) = E(7r/2,0) + [E aa (7r/2,0)(Sa) 2 +E /sp (7r/2,0)(S/3) 2 ]/2. (5) 

Note that the first derivatives E a and Ep vanish at a = 7tl2, f3 = since it is an easy 

direction, and the coordinate system has been chosen properly so that E a p{ftl2,0i) = 0. 

From (5) it is clear that E(a,P) possesses a high degree of symmetry and is completely 
determined by EaJjiflfS) and Epp(7c/2,Q>). This is consistent with the properties of Xin- 
A 3-d plot of Xin with E aa (7i/2,0) I Epp{it/2,Qi) = 1/2 is shown in Fig. 2 (where radii 
represent the values of /;„). 

III. Initial switching field H sw q along an easy direction and the 
correlation between H sw q and / in 



We now calculate the initial switching field of the particle along the x direction. This is 
done by requiring F aa (ft/2,0)F /3/3 (ft/2,0) - F a/} 2 1 2,0) = and solving this equation 
for H. Using (1) and E aj3 (ft/ 2,0) = , we obtain the result 

H swQ _ = [E aa (K 1 2,0) + E„ (7t 1 2,0)- I E aa {k 1 2,0) + E„ {k 1 2,0) I] / 2M s 

= min{E aa (x/2,0),E /j/3 (x/2,0)}/M s , (6) 
where rmn{E aa (ft/2,0),E /j p(ft'/2,0)} means the minimum of E aa (7ll2,0) and 
E^/2,0). 

To find the correlation between H sw q~ and the initial susceptibility we calculate the 
maximal value of From (4), Xin attains the maximum when e 1 u y = 1 or e 1 m = 1 , 
depending on the magnitude of E aa {7tl2,0) and E ^ (kI 2,0) . The maximum of Xin is 

Z„,max = M / / min{£ aa (^/ 2,0), ^ (^/ 2,0) } . e 2 Hy =1 or e l m =\ means that the 

external field H is either along the v-axis or the z-axis, which is perpendicular to the x- 
axis. This is consistent with the fact that M initially pointed to the positive x direction, 
since H perpendicular to the jc-axis gives maximal torque. From (6) and the definition of 
the intrinsic pinning field H p0 = M s I x in max we immediately obtain H p0 = H swQ _ . 



The key to understand the complementary nature of coercivity enhancement and 
exchange bias in the F-AF exchange coupled systems is the coherence of magnetization. 



The proof has been based on the coherent rotation of the F magnetization. Thus the 
complementary nature of coercivity enhancement and exchange bias confirms the 
rotation nature of the F magnetization in the F-AF systems. In fact, from the proof, the 
relationship H = M s I jtr inmax always holds for a single domain particle. In that sense, 

the complementary nature of coercivity enhancement and exchange bias in the F-AF 
exchange coupled systems is a special case for this relation. H sw q is the minimal reverse 
field to cause the magnetization to rotate. On the other hand, Xin,max corresponds to the 
maximum rotation of magnetization under an (infinitesimally small) external field 
applied at some direction perpendicular to the magnetization. The external filed is not 
applied in the same direction in these two cases. However, the relation 
H swQ - = M s I Xin max implies the same physical process. In other words, for a reverse field 
H = H the magnetization rotates in the direction that would attain the maximal initial 

swO ° 

susceptibility. When the external field is applied opposite to the magnetization, Zeeman 
energy has axial symmetry with respect to the magnetization direction. At H = H swQ _ , 

the magnetization is no longer stable and rotates in the direction that requires minimal 
energy increase. On the other hand, when the rotation of magnetization costs minimal 
energy, it attains maximal initial susceptibility. This is evident from 

Z„, m ax =M s 2 /mm{E aa (x/2,0),E^(x/2,0)} and (5). 



IV. The difference between the initial switching field and the actual 
switching field 

We distinguish between the initial switching field H sw q and the actual switching field. 
The initial switching field is defined as the reverse field at which the magnetization 
becomes unstable. In F-AF systems, sometimes it is found experimentally that the initial 
switching field is not the actual switching field for the F magnetization. That is to say, 
the F magnetization does not switch yet at H = H sM _ . Instead, it switches at a slightly 

larger field H sw . Fig. 3(a) shows one such example. In the graph, H sw0 ~ is the field at 
which the hysteresis loop begins to "curve" from the horizontal line of positive 
saturation. Note the actual switching field, at which the F magnetization changes to the 
reverse direction abruptly, is larger than the initial switching field. Appendix B shows 
that this happens when the overall anisotropy energy is asymmetric with respect to the 
stable direction of the magnetization in the absence of external field. To study the 
switching behavior, we need to know the energy profile for H > H swQ _ and we take the 

layered F-AF system as an example. For simplicity, we assume that the overall 
anisotropy of the F layer includes only a unidirectional component and a uniaxial 
component. Other cases can be discussed similarly. The easy direction of the 
unidirectional component of the overall anisotropy (denoted by e.d., also called the 
pinned direction) is assumed to differ from the easy axis of the uniaxial component 
(denoted by e.a.), shown in Fig. 3(b). Suppose without any external field, the F 
magnetization initially stabilizes at some intermediate direction between e.d. and e.a., 
corresponding to an energy minimum state. In the figure, the external field is applied 



opposite to the stable direction. Thus the overall anisotropy of the F layer is not 
symmetric with respect to that stable direction and the energy barriers have different 
heights on the left and right hand sides of the F magnetization (see Fig. 3(c)). If a field H 
opposite to the stable direction with increasing magnitude is applied, the location of the 
minimum in the energy profile does not change for H < H . But for H > H , the 

bJ r b swO swO ' 

above direction is no longer an energy minimum state. However, due to the asymmetry 
of the energy, a local minimum exists near that direction. Or equivalently the energy 
minimum shifts away from that direction continuously. This shift becomes larger as H 
increases. In this particular example, the shift is always away from the e.d. direction. 
Consequently, the magnetization rotates away from that direction continuously. No 
switching occurs until at some stage when H is so large that no energy minimum exists 
near that direction. When this happens, the magnetization will switch and the 
corresponding field is the actual switching field H sw ~, which is greater than H sw q. At 
H = H the F magnetization is unstable. More precisely, the state is metastable. It 

corresponds to a kink point in the energy profile. It is clear that the asymmetry of the 
energy profile has a great impact on the magnetization behavior. This is quite different 
from the symmetric case, where the pinned direction is the stable direction under 
H < H . and for H > H „ the pinned direction is unstable and corresponds to a 

swO ' .5vvO r r 

maximal point in the energy profile. No energy minimum exists near the pinned 
direction for// > H . In that case // SM T coincides with H sw q, resulting in a square 

hysteresis loop. In the asymmetric case a loop curved between H sw ~ and H sw q is 
observed, as shown in Fig. 3(a). 



The above argument can be verified by considering the free energy of the F-AF layered 
system. To be consistent with the notations in the previous sections, let Si and S 2 be the 
angles of e.d. and e.a. and a be the angle of the F magnetization, all from the stable 
direction of F magnetization in the absence of external field. When an external field is 
applied opposite to the stable direction, the free energy density of the system is 

F(a) = M S H cos a- K ed cos(« - S l ) + K ea sm 2 {a-S 2 ) 
where K e d and K ea denote the anisotropy constants for the unidirectional component and 
the uniaxial component, respectively. The derivatives of F(a) are 

F\a) = -M S H sin a + K ed sin(a -5 X ) + K ea sin 2(a -S 2 ), 

F" (a) = -M s H cos a + K ed cos(a - S i ) + 2K ea cos 2(a -S 2 ), 
F"\a) = M s H sin a - K ed sin(a - 8 X ) - 4K ea sin 2(a -S 2 ), 
F (4) (a) = M S H cos a - K ed cos(a -5,)- SK ea cos 2(a -8 2 ). 
a = is the stable state for H < H . At H = H the first and second derivatives of 

swO swO 

F(a) at a - vanish and the magnetization becomes unstable. Thus we have 
F'(0) = -K ed sin S i - K ea sin 28 2 = , 

F"(0) = ~M s H swQ _ + K ed cos S, + 2K ea cos 2 8 2 = . 

To study the energy profile near a - for H — H , we also need to calculate these 

higher order derivatives of F(a): 

F"'(0) = K ed sin S l + AK ea sin 28 2 



= -[-K ed sin S x - K ea sin 2 8 2 ] + 3K ea sin 28 2 
= -F'(0) + 3K ea sm 28 2 
= 3K ea sin 2S 2 , 
F (4) (0) = M s H m0 _ - K ed cos 8 X - SK ea cos 28 2 

= -[-M s H swQ _ + K ed cos 8, + 2K ea cos 28 2 ] - 6K ea cos 28 2 
= -F"(0)-6K ea cos2^ 2 
= -6K ea cos 2^ 2 . 

If 8 2 ^ we have F"'(0)^0 and a = is a kink point of the free energy. 
Consequently, an energy minimum near a - exists for H swQ _ < H < H . On the 

other hand, if 8 2 =0 we have F m (0) = and F (4) (0) = -6K ea <0. Thus a = is a 
maximum point of the free energy for H > H 



Fig. 3(c) shows schematically the changes of energy profile for the asymmetric case, 
especially for the reverse field near H sw q~ and H sw ~. In the figure, a is the angle from the 
stable direction. Note the shift of the location of the energy minimum and the change of 
energy profile near H sw0 ~ and H sw ~. It must be pointed out that both for symmetric and 
asymmetric cases, the relationship H = H holds strictly, as proved in section III. 



V. Conclusion 



Complementary nature of coercivity enhancement and exchange bias is generalized from 
the layered ferro-antiferromagnet (F-AF) exchange coupled systems to F-AF systems 
with arbitrary configuration. It is proved based on the coherent rotation of the F 
magnetization. By absorbing the effect of F-AF coupling into the anisotropy of the F 
part, we are able to avoid the complicated nature of the F-AF coupling at the interface 
and keep the generality of the discussion. A general discussion for the initial 
susceptibility of a single domain particle is provided. The fundamental correlation 
between the maximal initial susceptibility and the initial switching field along an 
arbitrary easy direction of a single domain particle is analyzed. The difference between 
the initial switching field and the actual switching field along an arbitrary easy direction 
is also discussed. 
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Appendix A 



Now we give a mathematical proof of the claim that Eo^tj/1,0) vanishes for some 
appropriate choice of the y and z axes when the positive x direction is an easy direction. 
Without loss of generality, suppose E a/} (ft/2,0) ^ in the current coordinate system. 

After rotating the y and z axes in the y-z plane with an angle y, we obtain a new 
coordinate system xy' z' (see Fig. 4). The y' and z' components of the magnetization are 

denoted by M , and M z , . Thus we have 
M y = cos jM , + sin ?M 7 , 
M z = -sin }M y +cos JM Z : . 

The v and z components of the magnetization M y and M z are expressed in terms of a and 
/3 by M = M s sin asin f5 and M , = M s cos a . From that we can obtain 

sin a sin j8 = cos ^sin a' sin j8'+ sin ^cos a' , 

cos « = - sin y sin a' sin cos y cos a' , 
where a' and ft are the polar and azimuthal angles of the magnetization in the new 
coordinate system. 

The above two relations give the dependences of a and /? on a' and ft . Taking partial 
derivatives with respect to a' and ft on both sides of these equations yields 



cosasin f3da Ida'+smacos /3df3lda'= cos ^cosa'sin /?'-sin ysma' , 
cosasm/5dal 3/?'+ sin a cos J3dfi ld/3'=cosy sin a' cos/3', 
-sinada/da'= -sin ^cosa'sin fi'-cos ys'ma' , 
- sin ad a I d/3' = - sin y sin a' cos /?' . 



Note after the rotation the x-axis remains unchanged. For the positive x direction, 
a = K ll,/3 = and hence a'= 7tl2,/3'= 0. Thus we have 

a^/aa'l (jr/20) = -sin7, (Al) 

tyi*p\ WVSi =fx*y, (A2) 

-3a/3a'l (;r/20) =-cos^ or 3a/3a'l (jr/20) = cos 7, (A3) 
-da/dj3\ /2Q) =-smy or aa/a^'l (jr/20) = sin 7. (A4) 
On the other hand 

E B . = dE I da' = dE Idadalda'+dE I dfidplda' = E a da/da'+E p d/3 /da', 

and 

E a ,p = d 2 Elda'dp'= (E aa da/djB'+E ap d/3/d/3')da/da'+E a d 2 a/da'djB' 
+ (E^da/djB'+E^djB/djB^djB/da'+E^jB/da'djB' . 
Since the ^-axis is an easy axis, E a {7l I '2,0) = 0, E p {7C I '2,0) = 0. Assuming E aj} = E pa 
and using (Al) - (A4), we obtain 

E a ,p, (7t 1 2,0) = [sin yE aa {k 1 2,0) + cos yE af3 (n 1 2,0)] cos y 

+ [sin yE afs (7tl 2,0) + cos yE pfs (71 1 2,0)] (-sin y) 



= sin ycos y[E aa {n 1 2,0) - E pp (Ttl 2,0)] 

+ (cos 2 y- sin 2 Y)E a/3 (K /2,0) 
= sin 2y[E aa (x/ 2,0) - (;z7 2,0)] + cos 2yE ap {Ttl 2,0) 
We choose y such that 

(* / 2,0) 2,0) 

cot 2r = — ^ . 

2^(^/2,0) 

Then E a7? , (;z7 2,0) = . 



Appendix B 

In this appendix, we derive a necessary and sufficient condition for H w ^ ^ swQ - f° r a 

single domain particle and discuss a method to determine H sw0 ~ from the hysteresis loop. 
The third partial derivatives of the overall anisotropy are assumed to exist in a 
neighborhood of the initial F magnetization direction. 



(1) Necessary and sufficient condition for H ^ H for a single domain 
particle 



H sw q has been defined as the reverse field along an easy axis at which the magnetization 
along the easy axis becomes unstable. By "unstable", it does not necessarily mean that 



the magnetization will switch immediately. The magnetization direction can still undergo 
smooth changes afterwards. Whether it switches or not depends on the energy profile. 
Thus we need to examine the energy profile and find the local minima for the system at 
H > H . from which the condition for H ^ H „ can be obtained. 

swi) sw swO 



As before, let a and /? be the polar and azimuthal angles for the magnetization. Suppose 
that the +x direction, i.e., a — nil and /? = , is an easy direction of the magnetization 
and H is applied along the -x direction. The magnetization stays in the +x direction when 
H < H swQ _ . For H > H , the magnetization becomes unstable. Our task is to find the 

stable magnetization direction when H > H „_ . For that, we increase an infinitesimal 

amount of the field 8ti > from H sw q and find the condition for a stable state to exist in 
a neighborhood of the +x direction. With H = H „ + SH we seek solutions to 

F a =dF/da = 0, (Bl) 
Fj8 = dF/dj8 = 0. (B2) 

Under what conditions can we find a stable solution in a neighborhood of the +x 

direction? 



We discuss a necessary condition for that first. Note F a =E a + M S H cos a cos J3 and 
F p =E p - M S H sin asin j3 . For I 8a \« 1 and I 5/3 \« 1 , 



F a {xl2 + 8a,Sft) = Y d ™ ^ 



v=0 



([Sa^- + S/3^ Q r +l F a )(x/2 + rjSa, rjSfi) 

+ ^ ^ , (B3) 

(n + l)\ 



where 77 e [0,1]. 



Suppose that the coordinate system has been selected properly so that 
E a p{7t 1 2,0) = Ep a {7t 1 2,0) = 0. Without loss of generality we assume 

E aa {71 1 2,0) < Epp (jl 1 2,0) . From our result in section III, 

H sw0 - =mm{E aa {7tl2,0\Epp{7tl2fi)}IM s = E^Tt 1 2,0)1 M s . (B4) 

From (Bl), (B3), and (B4), we obtain for SFL > , 

F a {nl2 + Sa,SP) = F a {7zl2,0) + {[Sa^ + Sp^AF a ){xl2,0) 

da dp 

([ Sa^- + Sj3^-] 2 F a )(x/ 2,0) 
da dp 

2! 

([Sa^ a + Sp^fF a )(x/2 + ri l Sa,ri 1 SP) 

+ _ 

= E a {7Tl 2,0) + [E aa {71 1 2,0) -M s H]Sa + E a p{7T/ 2,0)50 
+ ^[E aaa {Til 2,0){8a) 2 + 2E aa p {n 1 2,O){8a){80) 
+ E a pp{7r/2,0){S/3) 2 ] 



+ _ 

= E a {7tl 2,0) + [E aa {tzI 2,0) -M s (H swQ _ + M)]Sa 

+ E ap {7t 1 2,0)SJ3 + X -YE aaa {7C 1 2.0X&*) 2 

+ 2E aa/3 {7i 1 2,Q){da){dfi) + E m {nl 2,0)(#?) 2 ] 

{[Sa^- + S/3^-] 3 F a ){7u/2 + j 7l Sa,j] l SP) 
da dp 

3! 

= ^[E aaa {7T/2,0){Sa) 2 + 2E aap {k I ' 2,Q){Sa){dfi) 
+ E afsp {7tl2,Q){SP) 2 ] - M , {M){5a) 
(^SaJL + Sp^f F a ){7u/2 + Tj 1 Sa,rj 1 SP) 

dC( dp 

3! 

= 0, (B5) 
where r} 1 e [0,1]. In (B5), E a {7Cl 2,0) = since the +x direction is an easy direction. 
Similarly from (B2) we obtain 

Fp{x/2 + da, dp) = F p {ttI 2,0) + {[Sa^- + Sj3 ^-\F p )(* 1 2,0) 

aa dp 

+ ^{[Sa^ + Sp-^?F p ){7tl2+ri 2 Sa,ri 2 SI3) 
= F fj {7tl 2,0) + F p a {x 1 2,0)Sa + F {x 1 2,0)3(3 

+ ^{[Sa^ + Sp^\ 2 F fi ){xl2+ri 2 Sa,ri 2 Sp) 
= E p {x 1 2,0) + E pa {x 1 2,0)Sa + [E {x 1 2,0) -M s H\S(3 



+ \([Sa^- + Sj8^-] 2 F / ,)(7r/2 + j j2 Sa,j j2 Sj8) 

2 aCC dp 
= [E pp (*/ 2,0) - M s (H swQ _ + M)W 

+ ^([Saj- + S^^] 2 F p )^/2 + Tj 2 Sa,rj 2 S^ 

= [E pp (;r/2,0) - E aa {Ttl 2,0) -M s M]Sj8 



+ -{[da— + d/3—fF p ){7tl2 + r] 2 da,T] 2 dfi) 
2 da dp 

= (B6) 
where 7 2 € [0,1] . Again E p (ft/2,0) = since the direction is an easy direction. 



Depending on the magnitudes of E aa {7ll2SS) and E pp (ft/2,0) we consider two different 



cases: 



(i) E(x/2,0)<E BB (x/2,0) 



In this case we obtain from (B6) 

([Sa-?- + Sj3^-] 2 F p )(x/2 + j j2 Sa,T j2 S/3) 
2[E m (x/2,0)-E aa (x/2,0)-M s SH] 

Or equivalently 



SP„ F m {7cl2 + 7j 2 Sa,Tj 2 Sp) 

8a 2[E ^ 1 2,0) -E aa (7T 1 2,0) -M S 8H] P 

= [F R „„(7t 1 2 + 7i-8a,7i 8B)8a 

2[E pp (x 1 2,0) -E aa (x 1 2,0) -M s 8H] ^ 

+ 2F pap {7tl2 + r] 2 8a,r] 2 8P)8p\. 
Taking the limit as SH — > (for necessary condition, 8a — > , 8/3 — > ), we obtain 
8/3 1 8a = 0. 

Therefore in the limit as 8H — > , (B5) is simplified to 
E aaa (7tl2,Q)-2M s 8Hl8a = 0, 

or 

E aaa (x/2,0)8a/8H-2M s =0. 
In order for 8a 1 8H < °° , we need E aaa (Til 2,0) ^ , and we obtain 

8al8H = 2M s l E aaa (n 1 2,0) . (B7) 
Thus E aaa {7l 1 2,0) ^ is a necessary condition for (B5) to have solutions, i.e., a 
necessary condition for a stable sate to exist at the neighborhood of the +x direction. 



In fact, E aaa (K 1 2,0) ^ is also the sufficient condition. To see it, notice when 
SH ^0, 

F aa {7tl2 + 8a, 8/3) = E aa {tc 1 2,0) + E aaa {k 1 2,0)8a + E aap {jc 1 2,0)8/3 

-M s H+±-([Sa^- + S/3^-] 2 F aa )(x/2 + 7] 3 Sa,T] 3 S/3) 
2 da dp 

= E aa {k 1 2,0) + E aaa {n 1 2,0)8a + E aap {n 1 2,0)8/3 



-M s (H sw0 _ + SH) 

2 da dp 
= E aaa {n 1 2$)8a + E aa/3 {n 1 2,0)8/3 ~M S 8H 

~ M s 8H 

where 7] 3 e [0,1]. In the above derivation, (B4) and (B7) were used. 



Similarly we have 

F ap (7T/2 + 8a, Sfi) = E ap {tcI 2,0) + E afSa (x/2,0)5a + E app {n 1 2,0)8/3 

+ ^{[8a^ + 8/3^] 2 F aa ){7r/2 + Tj 4 8a,rj 4 8/3) 

= E apa {n 1 2,0)8a + E app {n 1 2,0) 8/3 

+ ^{[8a^ + 8/3-^] 2 F aa ){x/2 + T j4 8a,j j4 8/3) 

~E apa {7r/2,0)Sa, 

where rj 4 e [0,1] , and 

F pp {X/2 + da, 8/3) = E pp {tt/2,0) - M s H 

+ \{[Saj a + dp-^?F aa ){7ul2 + r l5 8a,7 l5 df3) 

= E pp {7tl2,0)-M s {H sw0 _+8H) 



= E pp {71 1 2,0) - E aa {7t 1 2,0) - M S M 

+ U[Sa^- + S J 8^-] 2 F aa ){7T/2 + r j5 Sa,j j5 S/3) 
2 da dp 

~ E pfS {7C 1 2,0) -E aa {7C 1 2,0), 

where 7] 5 e [0,1]. 

To check the stability of the state a-7t/2 + Sa, f3 = 8f5 , we consider the sign of 
F aa {7t 1 2 + 5a,d/3)F pi) {7[ 1 2 + 5a,d/3)-[F al) {7[ 1 2 + 5a,5P)f as M ^0. 

F aa {7T/2 + da, SjB)F pp {tt/2 + da, 5/3) - [F ap {tt/2 + da, 5/3)] 2 

~ [Ep p {7U 1 2,0) -E aa {7U 1 2,0)]M s M-[E a p a {7U 1 2,0)daf 

~ \E pp {7t 1 2,0) - E aa {7cl2,Q)]MM > 
since SH > and E aa {7t 1 2,0) < E pp {7t 1 2,0) from the assumptions. The second order 
terms were dropped. 

Therefore a = 7t I2 + Sa , /3 = 8/3 , is a stable direction for the magnetization at 
H = H wQ + 5H and we proved that E aaa {7C / 2,0) ^ is the necessary and sufficient 
condition for it. The case of E aa {7t 1 2,0) > E pp {7C 1 2,0) can be discussed similarly and in 



that case, E ^{k I '2,0) t is the necessary and sufficient condition for the existence of a 
stable direction for the magnetization at H = H sM _ + SH . 



E aaa (x/2,0) ^ means that the anisotropy energy is asymmetric with respect to the xz- 

plane and E^(^/ 2,0) ^ means that the anisotropy energy is asymmetric with respect 

to the xy-plane. In the case of layered F-AF systems, the above conditions reduce to the 
anisotropy energy asymmetric with respect to the stable direction of the F magnetization 
in the absence of external field. 



(ii) E(7tl2,Q) = E B A7tl2,Q) 



This case is more complicated. In this case the first order terms in (B6) vanish and we 
must keep all the second order terms. Hence, we need to solve two quadratic equations 
(B5) and (B6) jointly to obtain Sa/SH and d/3/SH. The existence and the nature of 
solutions depend on the third derivatives of anisotropy energy E (with respect to a and /?). 
Furthermore, even if a solution exists, whether or not it is a stable direction of the 
magnetization for H = H + 8H is still under questioning. We still need to use the 

criterion F aa {71 1 2 + da, SP)F^ (kI 2 + da, SjB) - [F aj3 (tt/2 + da, SjfJ)] 2 > to determine 
its stability, which again turns out to depend on the third derivatives of E (with respect to 



a and p). Thus, in general no simple specific conclusion can be drawn for this case. 
Since in reality, the condition E aa (ft/2,0) = (ft/ 2,0) is seldom satisfied, we do not 
provide further detailed discussion here. 

(2) Methodology to obtain H sw0 ~ when it is not the switching field 



H sw0 ~ accurately in experiments. We discuss the dependence of M x , the component of the 



applied in the -x direction). For a general 3-d F-AF system, the physical processes 
involved in this discussion are assumed to be reversible. For a 2-d F-AF system, this 
condition is not required. 

Note when AH is small, the magnetization direction will be close to the +x direction. 
M x = M s sinacos/? 

= M S sin0z72 + Aa)cos(A/?) =M S cos(Aa)cos(A/?) 



This question has significance in practice when H ^ H 



and we want to obtain 




in the +x direction, on H = H + AH when AH is small (recall that H is 




(B8) 



Assume E aa (7t 1 2,0) < E ^{ft 1 2,0) . From the previous section we have A/31 Aa~0. 
Thus from F a (7tl2 + Aafi) = we get 

E a {nl2 + Aafi) + M s H cos(x/ 2 + Aa) 

1 



~E a {Kl 2,0) + E aa (7C 1 2,0) Aa + - E aaa {n 1 2,0)( Aa) 



2 



-M s (H swQ _ + AH)sm(Aa) 
« ^(^/2,0)Aa + |£_(^/2,0)(Aa) 2 - M S (H sw0 _ + AH)(Aa) 

= ^E aaa (x/2,0)(Aa) 2 -M s (AH)(Aa) 



= 0. 

which yields 



Aa~ 2M - AH. 
E„(/!72,0) 



Using (B8) we obtain 



AM = M -M « — M^AaY « * -(AH) 

2 s [£_(;r/2,0)] 2 

Thus H SW Q can be obtained by using the above relation for the part of H > H within a 

small range (see Fig. 3(a)). 



FIG. 1. External field H and magnetization M in spherical coordinates. M initially 
pointed to the positive x direction, an easy direction of M. 

FIG. 2. Plot of initial susceptibility of a single domain particle with 

E aa (n/2,0)/E /vl (jr/2,0) = l/2,ia the unit of M s 2 1 E pp {7t I '2,0) . The radii represent 

the initial susceptibility values in different directions. In the plot, the coordinates 
for the x axis and the y axis are from -1 to 1 and the coordinates for the z axis are 
from -2 to 2. The magnetization initially pointed to the positive jc direction, an easy 
direction of the magnetization. 

FIG. 3(a). Hysteresis loop along the F magnetization stable direction in the absence 
of external field, of a Ni8oFe 2 o/FeMn F-AF coupled bilayer, measured by MOKE, 
(b). Schematic view of the easy direction of the unidirectional component e.d. 
(pinned direction), the easy axis of the uniaxial component e.a. of the overall 
anisotropy of the F layer, and the stable direction of F for an F-AF layered system in 
the absence of external field, (c) Schematic view of energy profiles for an F-AF 
layered system with asymmetry in the overall anisotropy, especially for the external 
field near H sw q and H sw . 



FIG. 4. Rotation of the coordinate system about the jc-axis by an angle y giving a 
new coordinate system xy z 
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